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Coherent acceleration of material wavepackets in modulated optical fields 
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We study the quantum dynamics of a material wavepacket bouncing off a modulated atomic mirror 
in the presence of a gravitational held. We find the occurrence of coherent accelerated dynamics for 
atoms beyond the familiar regime of dynamical localization. The acceleration takes place for certain 
initial phase space data and within specific windows of modulation strengths. The realization of the 
proposed acceleration scheme is within the range of present day experimental possibilities. 
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Enrico Fermi, in his seminal paper 'On the origin of 
cosmic rays', conjectured that "cosmic rays are originated 
and accelerated primarily in the interstellar space of the 
galaxy by collisions against moving magnetic fields" Q. 
Following this work, the possibility of accelerating parti- 
cles bouncing off oscillating surfaces Q was extensively 
studied in accelerator physics, leading to the develop- 
ment of two major models: The Fermi-Ulam accelerator, 
which deals with the bouncing of a particle off an oscillat- 
ing surface in the presence of another fixed surface par- 
allel to it; and the Fermi- Pustyl'nikov accelerator, where 
the particle bounces off an oscillating surface in the pres- 
ence ofgravity. In the case of the Fermi-Ulam accelera- 
tor IS 3 it was shown that the energy of the particle re- 
mains bounded and the unlimited acceleration proposed 
by Fermi is absent Q. In the Fermi- Pustyl'nikov acceler- 
ator, by contrast, there exists a set of initial data within 
specific domains of phase space that result in trajectories 
speeding up to infinity. 

In recent years, the efficient transfer of large momenta 
to laser-cooled atoms has become an important problem 
for a number of applications such as atom interferometry 
and the development of matter-wave based inertial sen- 
sors in quantum metrology. Possible schemes of matter- 
wave acceleration have been proposed and studied. For 
example, a spatially periodic optical potential applied at 
discrete times to an atom 7] was found to accelerate it 
in the presence of a gravitational field 8] . This r5-kicked 
accelerator operates for certain sets of initial data that 
originate in stable islands of phase space. 

In this Letter, we propose an experimentally realisable 
new technique to accelerate a material wavepacket in a 
coherent fashion. It consists of an atom optics version 
of the Fermi- Pustyl'nikov accelerator 9] where a cloud 
of ultracold atoms falling in a gravitational field bounces 
off a spatially modulated atomic mirror. This scheme 
is different from previous accelerator schemes in the fol- 
lowing ways: i) The regions of phase space that support 
acceleration are located in the mixed phase space rather 
than in the islands of stability (or nonlinear resonances); 



ii) The acceleration of the wavepacket is coherent; Hi) It 
occurs only for certain windows of oscillation strengths. 

Our starting point is the analysis by Saif et al. 
that establishes the dynamical localization of atoms in 
the Fermi-Pustyl'nikov accelerator (or Fermi accelera- 
tor for short) and shows a diffusive behavior both in 
the classical and the quantum domains beyond the lo- 
calization regime [10]. We now extend these results to 
identify conditions leading to the coherent acceleration 
of the atoms. We find clear signatures of that behav- 
ior both in an ensemble of classical particles and for a 
quantum wavepacket. In quantum mechanics, however, 
the Heisenberg Uncertainty Principle restricts the phase- 
space size of the initial atomic wavepacket which may 
result in the coherent acceleration occurring on top of a 
diffusive background. 

We consider a cloud of laser-cooled atoms that move 
along the vertical indirection under the influence of grav- 
ity and bounce back off an atomic mirror This 
mirror is formed by a laser beam incident on a glass 
prism and undergoing total internal reflection, thereby 
creating an optical evanescent wave of intensity I(z) = 
Io ex.p(—2kz) and characteristic decay length k~ l outside 
of the prism. 

The laser intensity is modulated by an acousto-optic 
modulator as ^(J 



I(z, t) = Iq exp(— 2kz + e smut), 



(i) 



where uo is the frequency and e the amplitude of mod- 
ulation. The laser frequency is tuned far from any 
atomic transition, so that there is no significant upper- 
state atomic population. The excited atomic level(s) 
can then be adiabatically eliminated, and the atoms be- 
have for all practical purposes as scalar particles of mass 
m whose center-of-mass motion is governed by the one- 
dimensional Hamiltonian 
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where p is the atomic momentum along z and g is the 
acceleration of gravity. 

We proceed by introducing the dimensionless posi- 
tion and momentum coordinates z = zuo 2 / g and p = 



2 




FIG. 1: (a) Varience of the momentum distribution, Ap, as a function of the modulation strength A, for an ensemble of 
particles initially in a narrowly peaked Gaussian distribution originating from the area of phase space that supports accelerated 
trajectories. The variance Ap (solid line) and average number of bounces (n) (dotted gray line) are displayed after an evolution 
time t = 300, for an ensemble of 10000 particles. The initial distributions in momentum and in coordinate space are gaussians 
centered at z — and p = 2tt 2 with Ap(0) = Az(0) = 0.1. The momentum distribution at t = 300 is illustrated in (b) 
for A = 37r/2. The broad, barely visible background results from the tails of the initial gaussian distribution outside the 
area of phase space that supports accelerated trajectories. The numerical calculations correspond to Cesium atoms of mass 
m = 2.2 x 10 _25 Kg bouncing off an atomic mirror with an intensity modulation of e = 0.55. The modulation frequencies extend 
to the megahertz range, and k _1 = 0.55 /im. 



poo/ (mg), the scaled time t = ut, the dimensionless inten- 
sity Vo = hw 2 VL e ft/ '(4mg 2 ), the "steepness" n = 2kg /uo 2 , 
and the modulation strength A = uo 2 e/(2kg) of the 
evanescent wave field, in terms of which the Hamiltonian 
takes the dimensionless form 

v 2 

H= (co 2 /mg 2 ) H = ^- + z + exp - Asint)]. (3) 

When extended to an ensemble of non-interacting par- 
ticles, the classical dynamics obeys the condition of in- 
compressibility of the flow , and the phase space distri- 
bution function P(z,p,i) satisfies the Liouville equation 
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dt 
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(4) 



where p = — 1 + kVq exp [— k(z — Asint)] is the force on 
a classical particle. 

In the absence of mirror modulation, the atomic dy- 
namics is integrable. For very weak modulations the in- 
commensurate motion almost follows the integrable evo- 
lution and remains rigorously stable, as prescribed by the 
KAM theorem. As the modulation increases, though, the 
classical system becomes chaotic. 

In the quantum regime, the atomic evolution is deter- 
mined by the corresponding Schrodinger equation 



p^_ 

2 



Vo exp [—k(z — A sin t)] 



t/> (5) 



where k = hco 3 / \mg 2 ) is the dimensionless Planck con- 
stant, introduced consistently with the commutation re- 
lation [z,p] = i(uo 3 /mg 2 )h = ik for the dimensionless 



variables z and p. We use Eqs. (@) and © to study the 
classical and quantum mechanical evolution of an ensem- 
ble of atoms in the Fermi accelerator. 

For very short decay lengths ft -1 and atoms initially far 
from the mirror surface, we may approximate the optical 
potential by an infinite potential barrier at the position 
z = Asincjt. In that limit the atoms behave like falling 
particles bouncing off a hard oscillating surface 0, • 

Classical version of the problem p demonstrats the ex- 
istence of a set of initial conditions resulting in trajecto- 
ries that accelerate without bound. Specifically, the clas- 
sical evolution of the Fermi accelerator displays the onset 
of global diffusion above a critical modulation strength 
\\ = 0.24 Q, while the quantum evolution remains local- 
ized until a larger value A M of the modulation |?J 0, ^( . 
Above that point both the classical and the quantum dy- 
namics are diffusive. However, for specific sets of initial 
conditions that lie within phase space disks of radius p, 
accelerating modes appear for values of the modulation 
strength A within the windows |(| 



sir < A < y/\ + (s7r) 2 , 



(6) 



where s can take integer and half-integer values for the 
sinusoidal modulation of the reflecting surface considered 
here. 

In order to compare the classical and quantum atomic 
dynamics within these windows in the present situation 
we calculate the width of the momentum distribution, 



Ap = y (p 2 ) — (p) , as a function of the modulation 
strength A > A n . In the classical case, we consider an 
ensemble of particles with a gaussian initial phase-space 
distribution P(z,p, 0) centered in a region of phase space 
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FIG. 2: Mirror images of the classical and quantum mechan- 
ical momentum distributions, P(p), plotted for (a) A = 1.7 
and (b) A = 2.4, after a propagation time t = 500. The spikes 
in the momentum distribution for A = 1.7 are a signature 
of coherent accelerated dynamics. The initial width of the 
momentum distribution is Ap = 0.5. 

that supports unbounded acceleration, and record the 
dispersion in momentum of that ensemble after a fixed 
propagation time. The corresponding quantum problem 
is treated by directly solving the Schrodinger equation for 
an initial gaussian wavepacket of zero mean momentum. 

Classically, the momentum space variance Ap remains 
small and almost constant for very small modulation 
strengths, which indicates the absence of diffusive dy- 
namics, but starts to increase linearly as a function 
of the modulation strength for larger values of A. As A 
is further increased, though, we find that the diffusion of 
the ensemble is sharply reduced for modulation strengths 
within the acceleration windows of Eq. (0), see Fig. [I] (a). 
Following each window the momentum dispersion grows 
again approximately linearly with A. We interpret the 
sharply reduced value of the dispersion as resulting from 
the non-dispersive, coherent acceleration of the atomic 
sample above the atomic mirror which is illustrated in 
Fig. [I] (b) and absent otherwise. 

In the quantum case the Heisenberg Uncertainty Prin- 
ciple imposes a limit on the smallest size of the initial 
wavepacket. Thus in order to form an initial wavepacket 
that resides entirely within regions of phase space leading 
to a coherent and dispersionless acceleration, appropri- 
ate value of the effective Planck constant is to be eval- 
uated, for example, by controlling the frequency uj [TEj . 
For broad wavepacket s, the situation that we consider 
in this Letter, the coherent acceleration manifests itself 
instead, both in the classical and the quantum cases, as 
regular spikes in the marginal probability distributions 
P(Pi~k) — J dxP(x,p,t) and P(x,t) = J dpP(x,p,t), 
which are absent otherwise. 

Figure illustrates the marginal probability distribu- 
tion P(p,t) in momentum space for (a) A = 1.7 and 
(b) A = 2.4, both in the classical and the quantum do- 
mains. In this example, the initial area of the particle 
phase-space distribution is taken to be large compared 
to the size of the phase-space regions leading to purely 
unbounded dispersionless acceleration. The sharp spikes 
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FIG. 3: Momentum variance Ap as a function of A for an 
atomic de Broglie wave (thick line) and an ensemble of par- 
ticles (dashed line) initially in a gaussian distribution after 
a scaled propagation time t = 500. The initial probabil- 
ity distributions, shown in inset, have variances Ap — 0.5 
and Az — 1. The gray line shows scaled diffusion coefficient 
D\/Dq as a function of A (arbitrary units). 




FIG. 4: (a) Square of the momentum variance (dark lines) and 
coordinate space variance (gray lines) as a function of time 
for A = 1.7. The coherent acceleration results in a breathing 
of the atomic wave packet, as evidenced by the out-of-phase 
oscillations of the variances, (b) Dynamics for A = 2.4, a 
modulation strength that does not result in coherent acceler- 
ation. Note the absence of breathing in that case. The two 
values, A = 1.7 and A = 2.4, correspond to the first maxima 
and first minima in Fig. Same parameters as in Fig. 



in time evolved momentum distribution, P(p, t) appear 
when the modulation strength satisfies the condition of 
Eq. (0, and gradually disappear as it exits these win- 
dows. These spikes are therefore a signature of the co- 
herent accelerated dynamics. In contrast, the portions of 
the initial probability distribution originating from the 
regions of the phase space that do not support acceler- 
ated dynamics undergo diffusive dynamics. 

We can gain some additional understanding of the dif- 
fusive behavior from the close mathematical analogy be- 
tween the system at hand and the kicked rotor model . 
It has been established mathematically ^| that for large 
modulation strengths, the diffusive behavior of classical 
systems described by the standard map displays a mod- 
ulated growing behavior. For large A, the diffusion coef- 



ficient is given by 

= D Q - J 2 (K) - Jl(K) + J 2 2 (K) + jf , 

„ (7) 

where K = 4A and L> = if 2 / 2 = 8A 2 [Ig and Ji, J 2 and 
J3 are first-, second- and third-order Bessel functions. 
Recent experiments by Kanem et al [T3 also reports such 
a behavior for the ^-kicked accelerator in the case of large 
modulations. 

A comparison between the classical behavior and the 
quantum momentum dispersion as a function of A is il- 
lustrated in Fig. 03 while the scaled diffusion coefficient 
D\/Do is shown in arbitrary units. It is interesting to 
note that the dispersio n exhibits maxima for oscillation 
strengths, A m = + (s7r) 2 )/2, which reside at the 

center of the acceleration windows 0, indicating that 
those trajectories that do not correspond to the phase 
space area supporting accelerated dynamics display max- 
imum dispersion instead. 

From the numerical results of Fig. [2 we conjecture that 
the spikes are well described by a sequence of gaussian 
distributions separated by a distance 7r, both in momen- 
tum space and coordinate space. We can therefore ex- 
press the complete time-evolved wavepacket composed of 
a series of sharply-peaked gaussian distributions super- 
posed to a broad background due to diffusive dynamics, 
such that 

00 

P(p) = ^ e -P 2 /^P 2 e -(f-"-) 2 /4^ ( 8 ) 

Tl— — OO 



where e << Ap, and AT is a normalization constant. 



Further insight in the quantum acceleration of the 
atomic wave packet is obtained by studying its temporal 
evolution. We find that within the window of accelera- 
tion the atomic wave packet displays a linear growth in 
the square of the momentum variance and in coordinate 
space variance as a function of time. Figure 0] illustrates 
that, for modulation strengths within the acceleration 
window, the growth in square of the momentum variance 
displays oscillations of increasing periodicity whereas the 
variance in coordinate space follows with a phase differ- 
ence of 180°. The out-of-phase oscillatory evolutions of 
Ap 2 and Az indicate a breathing of the wavepacket and 
are a signature of the coherence in accelerated dynamics 
as it disappears outside. As a final point we note that 
outside of the acceleration window the linear growth in 
the square of the momentum variance, a consequence of 
normal diffusion, translates into a t a law, with a < 1 
which is a consequence of anomalous diffusion. 
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